Correction: Genetic interactions among ADAMTS metalloproteases and basement membrane molecules in cell migration in Caenorhabditis elegans

[This corrects the article DOI: 10.1371/journal.pone.0240571.].


Introduction
Members of the ADAMTS family of secreted zinc metalloproteases have important roles in animal development. Most of these proteases degrade extracellular matrix components such as proteoglycans or collagens [1]. Nineteen ADAMTS genes have been identified in the human genome, and mutations in many of them result in hereditary diseases that are related to disorders of the extracellular matrix [2,3]. Multiple ADAMTS proteases often function in a common developmental process. For example, the functions of ADAMTS-5, -9 and -20 are required for interdigital web regression [4]. ADAMTS-9 and -20 are needed for closure of the palate and for craniofacial morphogenesis and neural tube closure through their function in ciliogenesis [4,5]. ADAMTS-5 and -15 act in myoblast fusion [6]. These enzymes appear to function in a partially overlapping manner. However, the roles shared by these ADAMTS proteases in development still remain elusive. Among five ADAMTS genes in C. elegans, gon-1 and mig-17 play essential roles in the development of the somatic gonad [7,8]. Both GON-1 and MIG-17 localize to the gonadal basement membrane (BM) [9,10]. In the mig-17 mutants, DTCs meander and stray, resulting in an abnormal gonadal shape. In contrast, in the gon-1 mutants, DTCs rarely migrate and the gonads remain small. We have isolated and analyzed the genetic suppressors for DTC migration defects in mig-17 mutants. Dominant gain-of-function mutations in the fbl-1/fibulin-1 and let-2/collagen IV α2 chain, which encode BM molecules, have been frequently isolated as suppressors [11,12]. The fbl-1(gf) mutations are amino acid substitutions in the second epidermal growth factor−like motif of FBL-1C, one of the two splicing isoforms, and FBL-1C protein is secreted by the intestine to be released into the gonadal BM in a MIG-17 activity-dependent manner [11]. Suppression of mig-17 by fbl-1(gf) mutations is dependent on the BM molecule NID-1/nidogen [12]. let-2(gf) mutations are associated with amino acid changes within the triple helix domain or the C-terminal non-collagenous domain (NC)1 of collagen IV. LET-2 protein is secreted from body wall muscle cells and DTCs and localized to the gonadal BM in a MIG-17 activity-independent manner. Unlike the case of fbl-1(gf), suppression of mig-17 by let-2(gf) mutations does not require NID-1 [12].
Genetic suppressor analysis of gon-1 mutants revealed that loss-of-function (deletion) mutations in fbl-1 can suppress the shortened-gonad phenotype of gon-1 mutants. Because fbl-1 deletion mutants also exhibit the shortened-gonad phenotype, GON-1 and FBL-1 act antagonistically to regulate gonad formation [13]. This genetic interaction is mediated by the control of collagen IV accumulation in the gonadal BM: GON-1 acts to reduce collagen IV levels, whereas FBL-1 acts to maintain collagen IV levels [14].
Although mig-17 and gon-1 mutants are phenotypically very different, they are both suppressed or enhanced by mutations in let-2 and fbl-1. In this study, we isolated novel suppressor mutations in emb-9 for mig-17 and in let-2 for gon-1 gonadal defects. Together with the previously isolated suppressors and enhancers, we investigated the consequences when suppressors for mig-17 were combined with gon-1 mutations, and when gon-1 suppressors and enhancers were combined with mig-17 mutations. We found that some of the emb-9, let-2 and fbl-1 mutations suppressed both mig-17 and gon-1, whereas others suppressed only mig-17 or gon-1. These results suggest that mig-17 and gon-1 have their specific functions as well as the functions commonly shared between them for gonad formation.

Microscopy
Gonad migration phenotypes were scored using a Nomarski microscope (Axioplan 2; Zeiss). Analysis of gonadal phenotypes was performed at the young-adult stage as described [20]. Although the gonadal phenotypes of some strains used are published, we reevaluated them in this study. The levels of EMB-9-mCherry were quantified as follows. For each sample, confocal images of a sagittal section of the DTCs were obtained with a Zeiss Imager M2 microscope equipped with a spinning-disk confocal scan head (CSU-X1; Yokogawa) and an ImageEM CCD camera (ImageEM; Hamamatsu Photonics). Using ImageJ software, we measured the fluorescence intensities along three drawn lines, each of which crossed the DTC BM; the average background intensities inside the gonad were subtracted from the peak values of the line scan, and the resulting corrected values were averaged.

Homology modeling of NC1 domains
Homology modeling was conducted based on the crystal structures of bovine collagen IV NC1 domains [24] using the SWISS-MODEL server. Ribbon diagrams were created and edited with the Swiss-Pdb Viewer software.

Isolation of novel suppressors of mig-17 and gon-1 mutants
The C. elegans hermaphrodite gonad arms extend to the anterior-right and posterior-left areas of the body cavity. The U shape of the gonad arms are generated by migration of the gonadal leader cells, the DTCs, over the body wall during larval development (Fig 1). gon-1 mutants exhibit shortened gonads due to the premature termination of DTC migration and are sterile. In contrast, mig-17 mutants show misshapen gonad arms due to the misdirected migration of DTCs, but they are still fertile. Although both of these genes encode ADAMTS family metalloproteases, the phenotypes shown by these mutants are very different.
To identify genes interacting with mig-17 and gon-1, we isolated novel suppressor mutations for DTC migration defects of mig-17 and gon-1 mutants using EMS mutagenesis. We isolated 11 suppressor mutants of mig-17(k174), a null allele that has a nonsense mutation in the pro-domain (S1 Fig). Two of these mutants, k204 and k207, were genetically mapped near the center of LG III and acted as dominant suppressors for mig-17 ( Fig 2B). Next-generation sequence analysis of k204 and k207 identified mutations in emb-9, which encodes the h1 subunit of collagen IV. Both mutations were amino acid substitutions in the C-terminal NC1 domain. We generated plasmids carrying these mutant alleles of emb-9 and introduced them into mig-17 mutants. The extrachromosomal arrays containing these plasmids partially rescued the gonadal phenotype of mig-17 mutants (Fig 2B), indicating that emb-9(k204) and emb-9(k207) mutations are causative for mig-17 suppression.
https://doi.org/10.1371/journal.pone.0240571.g002 (su1006) and sur-5::GFP are marker plasmids that result in the roller movement phenotype and GFP expression in all somatic nuclei, respectively. We isolated non-roller and GFP − fertile animals from the F 2 or F 3 generation of transgenic animals treated with EMS (S2 Fig). One of the two gon-1 suppressors, tk101, acted as a dominant suppressor of gon-1(e1254) and was genetically mapped to the right end of LG X. tk101was a recessive suppressor for the gon-1 (q518) null allele (Fig 2C). Next-generation sequence analysis revealed an amino acid substitution in the N-terminal region of the triple helical domain of let-2, which encodes the h2 subunit of collagen IV. The extrachromosomal array containing this let-2 mutant plasmid partially rescued the gonadal phenotype of gon-1(e1254) (Fig 2C), indicating that let-2(tk101) is the causative mutation for gon-1 suppression.
To understand the roles shared by these two proteases, we examined how the suppressors of mig-17 affect gon-1 and how the suppressors and enhancers of gon-1 affect mig-17 gonadal defects. We examined all the combinations of double mutants. The representative phenotypes exhibited by these double mutants are shown in Fig 4, and their phenotypic penetrance scored at the young adult stage is shown in Figs 5 and 6. In mig-17(k174) animals and in double mutants containing mig-17(k174), suppression was assessed by whether the U-shaped gonad was formed as in the wild type. In gon-1(e1254) animals and in the double mutants containing gon-1(e1254), suppression was assessed by whether the gonad arms reached the dorsal muscle.
Three mutations found in EMB-9 (k204, k207 and tk75) and one in LET-2 (k193) were localized to the C-terminal NC1 domain of these collagen IV molecules. Using SWISS-MO-DEL, we deduced the three-dimensional structures of the NC1 domains of EMB-9 and LET-2 based on the crystal structures of bovine collagen IV NC1 domains [24] (Fig 7). We found that three amino acid substitutions in EMB-9, which are separated from one another in the primary  Both fbl-1(tk45) and emb-9 (g23cg46) deletions are potential null alleles, as they are expected to introduce termination codons shortly after the sequence, were closely apposed in the three-dimensional structure. In the type IV collagen meshwork, triple-helical collagen molecules connect to one another through NC1-NC1 domain interactions. EMB-9(k207 and tk75) mutations were localized to the NC1-NC1 interface regions, and EMB-9(k204) was close to these interface regions, suggesting that these amino acid substitutions may affect physical interactions between two NC1 trimers. However, because these mutants were able to proliferate as homozygotes, their mutant collagen molecules were most likely successfully assembled into the functional network in the BM. It is interesting that the amino acid substitutions EMB-9(k204 and k207), which strongly suppressed gon-1 and mig-17, and EMB-9(tk75), which strongly enhanced gon-1 and mig-17, were closely localized in the three-dimensional structure.

Collagen IV accumulation in the DTC BM
Based on an immunohistochemical analysis, we previously reported that the reduced accumulation of collagen IV in the gonadal BM in animals with fbl-1(tk45), a null mutation, can be compensated by gon-1(e1254) [14]. To examine the amount of collagen IV accumulation in the BM quantitatively, we used a functional emb-9::mCherry fusion reporter [18]. Third larval stage animals in which their DTCs were at or shortly beyond the first turn were selected, and the fluorescence intensity of the BM surrounding the DTCs was measured (Fig 8A-8C). Among the six single mutants examined, we observed that the intensity was slightly higher in fbl-1(tk45), significantly lower in let-2(k196) and significantly higher in gon-1(e1254) as compared with wild type. When combined with fbl-1(k201 and tk45) or let-2(k196 and g25) mutations, the level of EMB-9-mCherry accumulation in mig-17 mutants was not affected except for the case of let-2(k196), in which a slightly lower accumulation was observed (Fig 8B). Because the three mutations fbl-1(k201), let-2(k196) and let-2(g25) suppressed mig-17 but fbl-1 (tk45) did not, the levels of EMB-9-mCherry accumulation were not correlated with mig-17 suppression. We then combined fbl-1(k201 and tk45) or let-2(k196 and g25) mutations with gon-1. We observed that the levels of EMB-9-mCherry accumulation in fbl-1(k201), fbl-1(tk45) and let-2(k196), all of which suppress gon-1, were significantly lower than that of the gon-1 single mutants, whereas the level was not affected in let-2(g25), which enhances gon-1 (Fig 8C). Thus, it is possible that the reduction in EMB-9-mCherry accumulation is indicative of gon-1 suppression.

Discussion
In the present study, we isolated novel suppressor mutations of gonadal defects related to mig-17 and gon-1 mutants. We identified alleles of emb-9 as mig-17 suppressors and an allele of let-2 as a gon-1 suppressor for the first time. We found that some of the emb-9 (collagen IV a1), let-2 (collagen IV a2) and fbl-1 (fibulin-1) mutations suppressed both gon-1 and mig-17, whereas others suppressed only gon-1 or mig-17. These results suggest that gon-1 and mig-17 have their specific functions as well as functions in common that relate to gonad formation. Probably, the loss of the gene-specific functions is the cause of the very different phenotypes of the gon-1 and mig-17 mutants. deleted region [11,17]. C. Summary of effects of fbl-1, emb-9 and let-2 alleles on mig-17(k174) and gon-1(e1254) mutants. Blue and red boxes represent suppression and enhancement, respectively. Numeric values representing strength of suppression or enhancement were calculated as follows. Percentages of anterior and posterior gonadal defects were averaged for each strain. Strength of suppression indicates percent recovery of the gonadal defect in double mutants relative to that of the mig-17 or gon-1 single mutants. Strength of enhancement indicates percent decrease of the ratio of non-defective gonads in double mutants relative to that of the mig-17 or gon-1 single mutants.
https://doi.org/10.1371/journal.pone.0240571.g003 We presented a diagram of expected mechanisms of the action for MIG-17 and GON-1 which will be discussed here (Fig 9). The fbl-1(tk45) null mutation suppressed gon-1 but enhanced mig-17. The gon-1 suppression is likely to be due to the reduction of collagen levels in the late larval stages that results from loss of fbl-1 activity [14], although the collagen levels in fbl-1(tk45) were not reduced in the mid-L3 stage when the DTCs make their first turn ( Fig   Fig 6. Percentage of gonad arms that failed to reach the dorsal muscle of gon-1(e1254) mutants in the presence of  fbl-1, emb-9 and let-2 alleles. Red and blue bars represent defects in anterior and posterior gonad arms, respectively. n = 60 for each experiment. P-values from Fisher's exact test comparing the double mutants with gon-1(e1254) animals: �� P < 0.01; � P < 0.05; NS, not significant.
The gain-of-function mutations of collagen IV emb-9(k204, k207) and let-2(tk101, k196 and k193) were potent suppressors of mig-17 and gon-1. We found that the BM collagen levels in let-2(k196) were significantly decreased and that let-2(k196) suppressed the dramatic increase in collagen accumulation in gon-1 animals. Although we did not examine the other gain-offunction collagen mutations, it is possible that the levels of BM collagen are similarly affected. In contrast, the gain-of-function mutation emb-9(tk75) enhanced both mig-17 and gon-1. We previously showed that the levels of BM collagen in animals expressing the mutant EMB-9 (tk75) h1 subunit can be maintained in the absence of FBL-1, which is otherwise required for the maintenance of BM collagen levels [14]. Therefore, it is likely that too much accumulation of collagen in the BM could be causative for both the mig-17 and gon-1 mutant phenotypes  9). Although we could not detect over-accumulation of collagen in the mig-17 mutant, this could be because of the insufficient sensitivity of our assay condition.
This model is not, however, consistent with the fact that the loss-of-function mutations of collagen IV, emb-9(g34 and g23cg46/+) and let-2(g25 and b246), which conceivably reduce the levels of BM collagen, all acted as enhancers of gon-1 even though they acted as suppressors of 9-mCherry. Arrowheads point to DTCs. Bar: 20 μm. The right panel illustrates the gonadal BM (brown). The fluorescence intensity of the DTC BM was quantified by averaging the peak values along three lines (green) that cross the DTC BM (see Materials and methods). B, C. Box-and-whisker plot of the fluorescence intensity of EMB-9-mCherry in the DTC BM in animals with wild-type and mutant fbl-1 and let-2 alleles and with those mutant alleles in combination with mig-17(k174) (B) or gon-1(e1254) (C); n = 20. Boxplots indicate the median and the interquartile range. Whiskers extend to the minimum and maximum values within 1.5 times the interquartile range. Points indicate outliers. P-values from Student's t-test are indicated: �� P < 0.01; � P < 0.05; NS, not significant.
https://doi.org/10.1371/journal.pone.0240571.g008  [8,25]. Although the levels of BM localized FBL-1 and NID-1 are low before the activation of MIG-17 (A), collagen degradation by MIG-17 promotes more FBL-1 accumulation (B, C). The BM localized FBL-1 then recruits NID-1, which facilitates directional DTC migration [12] (D). E. Model for GON-1 action. GON-1 likely localizes to the gonadal BM from the L1 stage because the null allele gon-1(q518) completely blocks DTC migration [26]. The GON-1 collagenase activity is inhibited by FBL-1 localized to the BM and controlled in appropriate levels. GON-1 also acts in the endoplasmic reticulum to transport integrin receptors from the endoplasmic reticulum to the Golgi, and therefore to the plasma membrane.
https://doi.org/10.1371/journal.pone.0240571.g009 mig-17. The loss-of-function emb-9(xd51) mutation also enhances gon-1 in the presynaptic bouton overgrowth phenotype [27]. In contrast, emb-9(g34) and emb-9(b189) can suppress the synaptic defect of gon-1 when emb-9; gon-1 double mutants are shifted up from 16 to 25˚C after completion of embryogenesis [28]. In this case, it is possible that collagen IV containing the temperature-sensitive mutant proteins EMB-9(g34) and EMB-9(b189) form intracellular aggregates and thus are not secreted [16], and therefore the BM accumulation of collagen could be considerably decreased. Our analysis using the EMB-9-mCherry reporter revealed that the collagen levels of let-2(g25); gon-1 or let-2(g25); mig-17 double mutants were not lowered relative to the respective gon-1 or mig-17 single mutant (Fig 8B and 8C). It might be possible that a subtle reduction in collagen accumulation is sufficient for amelioration of BM physiology in mig-17 mutants but instead worsens that in gon-1 mutants. The slight collagen reduction may promote NID-1 accumulation in the BM and suppress the mig-17 gonadal defect.
Why is gon-1 enhanced by collagen IV loss-of-function mutations? This seemingly contradictory phenomenon may be related to the dual function of GON-1. In addition to its predicted extracellular protease activity, GON-1 functions in the endoplasmic reticulum to transport secreted or membrane proteins from the endoplasmic reticulum to the Golgi. This transport activity is dependent on the C-terminal GON domain but is independent of protease activity [29]. For example, cell surface expression of the integrin receptors that are required for cell migration may be reduced and, therefore, the DTC migration activity in gon-1 mutants may be weakened (Fig 9E). The fact that RNAi knockdown of the integrin h subunit ina-1results in a short-arm gonad phenotype similar to that of gon-1 supports this possibility [30]. Because the remodeling of the BM is coupled with epithelial cell migration [31,32], reduced migration of DTCs could also lead to downregulation of BM remodeling, resulting in thick accumulation of collagen, which can further block DTC migration. Remodeling of the BM is likely to be mediated by the GON-1 protease activity.
We previously showed that a reduction in collagen IV lowers PAT-3/β-integrin expression in DTCs [14]. A similar correlation between collagen IV and the integrin expression level has been observed in pancreatic β cells migration on Matrigel [33]. At the early second larval stage, the fluorescence levels of EMB-9-mCherry in gon-1 mutants were closer to those in the wild type (K.N., unpublished). Thus, it might be possible that when combined with loss-of-function collagen IV mutants, the reduced collagen levels may further impair integrin expression in DTCs of gon-1 animals in the second larval or younger stages. This might be the reason for the enhancement of the gon-1 phenotype in the presence of collagen IV loss-of-function mutations.
It is unclear why the gain-of-function collagen IV mutation let-2(k196) did not enhance gon-1, as did let-2(g25), even though let-2(k196) also reduced collagen accumulation. We speculate that GON-1 is the major enzyme responsible for turnover of BM collagen IV and that the mutant LET(k196) protein may confer the property by which the collagen IV meshwork is turned over quickly as compared with the wild-type meshwork even with the weakened activity of GON-1(e1254) mutant enzyme. If this is the case, the collagen meshwork containing EMB-9(tk75) may have gained a slower turnover rate.
Our observations suggest that both MIG-17 and GON-1 function to reduce collagen IV in the BM. Although this is consistent with the idea that collagen IV is the direct substrate of these enzymes, we still do not have evidence of this interaction. Future biochemical analysis is thus needed to determine the substrates. In addition, molecular structural analysis of interactions among these ADAMTS proteases and the BM molecules should shed light on the detailed mechanism of BM remodeling during organogenesis.